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Abstract. It is widely thought that neutralinos, the Ughtest 
supersymmetric particles, could comprise most of the dark 
matter. If so, then dark halos will emit radio and gamma ray 
signals initiated by neutralino annihilation. A particularly 
promising place to look for these indicators is at the center 
of the local group dwarf spheroidal galaxy Draco, and recent 
measurements of the motion of its stars have revealed it to be 
an even better target for dark matter detection than previously 
thought. We compute limits on WIMP properties for various 
models of Draco's dark matter halo. We find that if the halo 
is nearly isothermal, as the new measurements indicate, then 
current gamma ray flux limits prohibit much of the neutralino 
parameter space. If Draco has a moderate magnetic field, 
then current radio limits can rule out more of it. These results 
are appreciably stronger than other current constraints, and 
so acquiring more detailed data on Draco's density profile 
becomes one of the most promising avenues for identifying 
dark matter. 



1 Introduction 

Despite the popularity of the neutralino as a dark matter can- 
didate, efforts to constrain its properties have only been able 
to rule out a small fraction of the relevant parameter space. 
Direct detection experiments based on scattering with nuclei, 
and indirect searches based on neutralino annihilation, have 
shared this difficulty. But recent results can improve the sit- 
uation. New observations of the Draco dwarf galaxy reveal 
that it is strongly dark matter dominated, and that its dark 
matter distribution is at least nearly isothermal. In this paper 
we investigate the detectability of WIMP annihilation sig- 
nals from several appropriate halo models. Interestingly, if 
Draco is indeed isothermal down to very small radii, we will 
show that current gamma ray measurements rule out a signif- 
icant fraction of the neutralino parameter space. If one makes 
an additional assumption that Draco harbors magnetic fields, 
then current radio measurements rule out additional param- 
eter space. For example, a 1 /iG field would eliminate most 
of it. Other halo models yield weaker current limits, some of 
which become significant with upcoming detectors. 

The choice of a particular dark matter distribution is crit- 
ical to both gamma ray and radio limits, and we devote sec- 
tion 2 of this paper to that consideration. The assumption 
of a magnetic field in Draco is discussed in section 4. The 
remainder of this introductory section provides the relevant 
background. 



The cold dark matter (CDM) scenario starts with a weakly 
interacting massive particle (WIMP), created abundantly in 
the big bang and surviving today to dominate the matter den- 
sity in galaxies. The velocities of stellar orbits imply that 
dark matter occupies galaxies with a particular radial density 
profile. Over a wide range in radius, the density appears to 
fall off as p '-^ r~^, which is the same structure as would be 
found in a self-gravitating isothermal sphere. But when sim- 
ulations of the cosmological CDM evolution are performed, 
taking into account the influence of A^-body dynamics and 
tidal interactions, the profiles are found to be more compli- 
cated at large and small r (e.g., Navarro, Frenk, and White 
(1996), Moore et al. (1999)). We discuss these profiles in 
more detail in subsection 2. 1 . 

Simulations universally generate a central cusp in every 
dark halo. Slopes vary, but p ^ remains for some 
positive 7, and therefore the density blows up at the center 
However, a controversy ensues at this point because obser- 
vations often differ. 21 cm rotation curves from low surface 
brightness (LSB) galaxies have presented evidence for cen- 
tral cores of approximately constant density (McGaugh and 
de Blok, 1998). A separate study of 16 LSB galaxies con- 
tends that beam smearing effects render these measurements 
unreliable, such that cusps are still fully consistent with the 
data (van den Bosch et al., 2000). Many of the same galax- 
ies have since been re-examined with high resolution opti- 
cal rotation curves, and that work found 30 LSB galaxies 
evidently possessing cores, not CDM cusps (de Blok, Mc- 
Gaugh, and Rubin, 2001). Dwarf galaxies too have long 
been considered problematic for CDM because they typi- 
cally have the linearly rising rotation curves of a constant 
density core (Navarro, Frenk, and White, 1996), but a subse- 
quent study of 20 dwarf galaxies gives rise to interpretations 
ranging from cores to steep r^^ cusps (van den Bosch and 
Swaters, 2001). These issues constitute an area of much de- 
bate, but in general, because many observations favor cores 
over cusps, it is of great interest to determine what processes 
might prevent or destroy a cusp. 

If a WIMP cusp were to realize, and if WIMPs exist to- 
day in equal numbers with antiWIMPs (or if they are their 
own antiparticles), then we can expect densities high enough 
for dark matter annihilation to become significant. If the an- 
nihilation products lead to observable particles, then cusps 
become hot spots to search for WIMPs. 

Perhaps the most likely WIMP candidate is the neutralino 
(x), which is the lightest supersymmetric particle (LSP) in 
the minimal supersymmetric extension of the standard model. 



Supersymmetry is a symmetry relating fermions to bosons, 
and thus it introduces a new class of particles which are "su- 
perpartners" to standard model particles. The theory helps 
to remedy the hierarchy problem set by the electroweak and 
Planck scales. In models with conservation of R-parity, a su- 
persymmetric particle is prohibited from decaying into only 
non-supersymmetric ones, and since the neutralino is the LSP, 
there is no supersymmetric state of lower energy available. 
Thus the x is stable, making it a well motivated dark mat- 
ter candidate (for a review of supersynmnetric dark matter, 
see Jungman, Kamionkowski, and Griest (1996)). 

Each XX annihilation produces quarks which bind primar- 
ily into pions. Neutral pions usually decay to gamma rays; 
charged pions usually decay to neutrinos and muons, and 
subsequently electrons (and positrons), which generate radio 
waves by synchrotron radiation in an ambient magnetic field. 
(This approach does not wholly rely on supersymmetry, but 
rather on some DM particle which annihilates into pions.) 
Therefore, a cusp produces observable radio and gamma ra- 
diation (Berezinsky, Gurevich, and Zybin, 1992; Berezinsky, 
Bottino, and Mignola, 1994); the trick is just finding a galaxy 
with a cusp. 

The candidate galaxy should have a very high mass-to- 
Ught ratio (M/ L) to ensure that its internal gravitation is DM 
dominated, and it should be nearby to facilitate the detection 
of a weak signal. As it turns out, the highest known M /L 
comes from a dwarf satellite of our own Galaxy. The Draco 
dwarf spheroidal galaxy, only 79 kpc distant from Earth (van 
den Bergh, 2000), was recently re-examined for detailed stel- 
lar velocities (Kleyna et al., 2001a). A new, higher than pre- 
viously estimated M/L was found to be 440 ± 24OM0/L0, 
with a mass of M ~ 8.6 x 10'' Mq. Draco has the mass of a 
dwarf galaxy and the starlight of a globular cluster. 

Furthermore, the candidate galaxy should have a high cen- 
tral WIMP density, without contaminants such as massive 
baryonic clouds or a central black hole. The stellar veloc- 
ity dispersion measurements in Draco are consistent with an 
isothermal dark matter halo (Kleyna et al., 2001a), implying 
a steep central density. Its present hydrogen gas population 
is less than 450 Mq (Young, 1999), with no evidence of a 
central black hole. The Sloan Digital Sky Survey reported 
that tidal deformation can only exist at a level of < lO^-^ of 
the central stellar surface density (Odenkirchen et al., 2001), 
supporting the case for dynamical relaxation and therefore 
meaningful extraction of mass density from velocity infor- 
mation. Given these properties, we present particle dark mat- 
ter constraints enabled by observations of the stars in Draco. 

The next section provides more detail regarding halos with 
cusps and cores. Section 3 follows with the calculation of 
the observable flux created by WIMP annihilations, includ- 
ing electron losses by inverse Compton scattering, which is 
an important correction not usually considered in such anal- 
yses. Section 4 presents the results of these calculations, us- 
ing observational bounds from EGRET and the VLA, with a 
discussion of anticipated improvements by new experiments, 
and a comparison among various halo profiles. We summa- 
rize in section 5. 



2 The Dark Halo 

A power law cusp's density profile has the form p ~ . 
Binney and Tremaine (1987) have put forth some reasoning 
for the expectation of a singular isothermal sphere (SIS) dis- 
tribution to a first approximation, where 7 = 2. But do steep 
cusps like this really exist in nature? 

Observations of galaxy centers with the Hubble Space Tele- 
scope (HST) shed some light on this question (Crane et al., 
1993; Gebhardt et al., 1996). These authors have constructed 
three dimensional luminosity densities based on an assumed 
spherical geometry in the target elhpticals and in the bulges 
of spirals, further assuming that the light traces the mass 
in these systems. Among power law galaxy centers, two 
populations emerge: bright galaxies (roughly My < —20) 
with slopes in the range 0.2 < 7 < 1.5, and faint galax- 
ies {Mv > —20) with 7 ranging between 1.5 < 7 < 2.5. 
The peaks in this bimodal distribution are at 0.8 and 1.9, re- 
spectively. Many faint galaxies in the sample therefore have 
nearly isothermal density profiles at their centers. 

2.1 Halo Models 

An ordinary kinetic-molecular gas, self-gravitating and held 
isothermal, will assume the SIS density profile p r~^. In- 
stead, we have a collisionless system of WIMP particles. If 
they are allowed to evolve under their collective gravitational 
potential, the initial infall can be expected to be clumpy, and 
as a result, the overall potential seen by any one WIMP is the 
superposition of the central field and the temporary group- 
ings. Trajectories are thus altered, and over a few orbital pe- 
riods, a new smooth centrally weighted distribution will be 
obtained. This process was termed "violent relaxation" by 
its author, Lynden-Bell (1967), wherein the central collapse 
is achieved by transporting angular momentum outward. 

The solution at the end of this relaxation turns out to be 
the Maxwellian velocity distribution with constant velocity 
dispersion at every radius, although the population dwindles 
somewhat beyond the virial radius, as those particles do not 
often interact with the lumpy center and therefore do not join 
the infall. One important caveat, however, is that the violent 
relaxation process proceeds until a steady state is achieved, 
whether or not that state is the Maxwellian result of complete 
relaxation. This means that it is possible for the final state 
to retain features unique to its own collapse sequence (one 
reason why simulations are useful for this work). 

Supposing the end state of the process is complete relax- 
ation, we have a system of self-gravitating particles with a 
Maxwellian velocity distribution at every point in space, the 
solution of which is again the SIS density profile (Binney 
and Tremaine, 1987). If the system fails to achieve com- 
plete relaxation, then the iimermost orbits may not be fully 
populated, and one might expect 7 < 2 as is often seen in 
numerical simulations. 

Rotation curves of galaxies usually show a constant ve- 
locity over a wide range of radial values. On the contrary 
however, some observations reveal core regions of about 1 — 
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10 kpc in radius which appear to enclose an r-independent 
matter density. This is evident in the rotation curves of spiral 
galaxies which tend to rise hnearly at small radii, and those 
of some dwarf galaxies which rise Unearly for the full extent 
of the observed galaxy. To model this, one typically invokes 
a modified isothermal profile: a spherically symmetric den- 
sity profile which smoothly blends a constant density core 
with isothermal behavior outside the core. 

In order to calculate observables from annihilation reac- 
tions, one must start by assuming some profile. For exam- 
ple, gamma rays from annihilating neutralinos were origi- 
nally calculated by Berezinsky, Gurevich, and Zybin (1992) 
using p ~ r~^-^ for the Galactic center. In Gondolo (1993), 
the same signature from the Large Magellanic Cloud was 
calculated assuming p ~ a^/(r^ + a^), a modified isother- 
mal model with core radius a. Non-constant central densi- 
ties with radial power law exponents have been employed 
in different contexts: 7=1 (NFW: Navarro, Frenk, and 
White (1996), used for example in Gondolo (2000) and Blasi, 
Ohnto, and Tyler (2002)), 7 = 1.5 (Moore: Moore et al. 
(1999), used for example in Calcaneo-Roldan and Moore 
(2001)), and 7 = 2 (isothermal: used for example in Blasi, 
Olinto, and Tyler (2002)). It is worth mentioning that 7 ap- 
pears to vary with the resolution of the CDM simulation, and 
at small radii, a reaUstic model of the high cusp density can 
be expected to require more accurate simulations. 

In this work, we consider seven distinct models of Draco. 
These include steep power-law cusps with 7 < 2, the Moore 
profile, and isothermal halos truncated by constant density 
cores. 

A power law DM halo profile must include a small con- 
stant density core inside the cusp for any 7 > 0, due to the 
fact that inside some radius, the annihilation rate gets so large 
that the over-density is destroyed as fast as new infall can 
fill the region. The constant density region, found by setting 
cusp forming time scale (the clump crossing time, or approx- 
imately the free fall time ~ 1 / \/Up) equal to the annihi- 



lation time scale (1/n^ 



{), establishes a characteristic 



inner radius (Berezinsky, Gurevich, and Zybin, 1992). For 
example, if 7 = 2 (SIS), we have (Blasi, OUnto, and Tyler, 
2002): 



«halo 



(1) 



where the xx annihilation rate (cross section x velocity) 
{(jv)-x^x is approximately constant in the low velocity Umit 
appropriate for galactic dark matter (Berezinsky, Gurevich, 
and Zybin, 1992), and the time scales mentioned above are 
~ 5 X 10^ years at -Rmin- The external radius iiext of Draco's 
halo is where it blends smoothly into the Milky Way halo. If 
we define the cuspy halo density as 
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(2) 



then i?ext and the coefficient A can be jointly set by equat- 
ing the X density at the edge of Draco's halo to the ambient 
Milky Way halo density at Draco's location, rihaio- Addition- 
ally, integrating eq. (2) out to i?ext must give Draco's total 



mass. The following are numerical values for these quanti- 
ties, taking = 100 GeV and (crw)xx = 10^^^ cm^/s, 
listed here for the 7 = 2 case (the generalization to 7 < 2 is 
straightforward): 



-Rext 
-^min 
?^halo 
n(-Rmin) 



4.6 kpc 



= 9.4 X 10" cm 
= 2.2 X 10"^ cm"^ 
= 5.7 X 10^ cm"^ . 



As will be discussed in subsection 3.1, the gamma ray flux 
observed from a DM clump depends on iJmin, so it is worth 
considering whether it is appropriate to use the free fall time 
as the cusp formation time. Since violent relaxation is a com- 
plex process one might imagine some longer characteristic 
time scale. Without any particularly compelling choice for 
this time scale, we can set a conservative upper Umit by using 
the age of the universe. For Draco, using Tuniv = 15 Gyr, the 
resulting gamma ray flux is found to be about 5 times weaker 
than it would be with eq. (1), for any choice of m-)^. 

The NFW and Moore halo profiles mentioned previously 
are the most popularly quoted "universal profiles" obtained 
from simulations. Both employ power-law cusps internal to 
some scale radius, outside which the density faUs off as r~^. 
The behavior is nearly isothermal in the vicinity of the scale 
radius . 

In calculating iiext. we have used the customary NFW 
profile for our own Galaxy, 



?^halo = "•0 



1 



r 



(3) 



whose scale radius ~ 30 kpc, total size i?ext — 300 kpc, 
and characteristic density no — 7.5 x 10"'' cm~^, can be 
set by the boundary conditions imposed by the halo density 
local to our solar system (6.5 x 10~^^ g/cm^), and the total 
mass of the Galaxy (10^^ Mq), subject to an estimate that 
Rexk/Ts — 10. As it turns out, various Galactic halo mod- 
els roughly agree on the local dark matter density at Draco's 
distance from the center, so the choice of NFW as the Milky 
Way DM profile is not particularly important here. 

The Moore profile (Moore et al., 1999) is more recent, de- 
riving from higher resolution computations than NFW. We 
therefore include it in our calculations of Draco. The Moore 
profile is 



n(r) 



no 



-1.5 



1.5 



(4) 



where no ~ 2.5 x 10^^ cm^'^ and r,, ~ 0.35 kpc can be fixed 
in similar fashion, using Draco's total mass, the Milky Way 
density local to Draco, and a concentration factor i?ext/''« — 
10 as boundary conditions. 

In a recent paper (Power et al., 2002), a convergence study 
is described which attempts to resolve central halo densi- 
ties in CDM simulations. The result for a galaxy-sized halo 
appears to be a decreasing exponent moving inward from 
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the virial radius (comparable to -Rext in this work), reach- 
ing 7 < 1-2, without any particular inner power law slope 
convergence. Because this is a state-of-the-art simulation, 
the fact that the inner slope is not especially steep should not 
be ignored. However, the authors note that their innermost 
resolved point is only 0.5% of the virial radius, which is out- 
side our region of interest for Draco, and that the effect of 
poor resolution is usually to generate artificially low central 
densities. The reasons which cause this simulation to differ 
in inner slope from others Uke Moore are currently unclear, 
but external to a few times the innermost resolved radius, the 
various CDM halo profiles are difficult to distinguish. In any 
case, the central drop in 7 is not evident in Draco. 

2.2 Modehng Draco 

The stellar motion measurements from Draco are best fit by 
7 ~ 1.7 over ~1 kpc (Kleyna et al., 2001a); strong veloc- 
ity anisotropies and strong deviations from isothermality are 
ruled out over this radial scale. (Star counts and luminosity 
density are insufficient for determining the mass profile, be- 
cause according to the same stellar velocty study, mass does 
not follow light in Draco.) In the central ^0.2 kpc, 7 appears 
to be > 2, although other cusp slopes and cores cannot be 
ruled out at small radii with the current data set. So although 
the observations support the SIS at intermediate radii, they 
don't directly speak for radii smaller than about 10 pc, inside 
which our annihilation signal originates. 

Since annihilation rate goes as the square of the density of 
X particles (n^ {<^v)^^ cm~^ s~^), the choice of 7 is very im- 
portant to the result. Observed 7 values from rotation curves 
are neither accurate enough nor consistent enough to use uni- 
formly for this purpose. Although a positive detection of 
gamma or radio signals could identify a particle dark mat- 
ter source, the difficulty in choosing 7 with any confidence 
has made it hard to realistically constrain WIMP parameters 
such as mass and cross section on the merit of observational 
upper limits alone. 

However, this situation may soon change. Draco's mass 
has recently been shown to be completely dominated by dark 
matter, with stellar velocities that favor nearly isothermal or- 
ganization, and the authors of these observations contend 
that the prospects are excellent for improving the velocity 
curve data further by sampling more stars (Kleyna et al., 
2001a). This may become a particularly fruitful approach 
to further particle DM studies. So to agree well with current 
data on Draco, we use halo models with 7 = 2.0, 1.9, 1.8, 
and 1 .7, the Moore profile, and isothermal profiles truncated 
with cores at 1 pc and 0.1 pc. 

Other galaxies currently have more detailed velocity dis- 
persion or rotation data. We choose Draco for this study 
partly for its proximity, and partly for its lack of features 
which would be expected to disrupt the cusp. Baryonic mat- 
ter and central black holes are factors which contaminate 
this kind of analysis in other galaxies including our own, but 
Draco is evidently little more than a spherical clump of self- 
gravitating dark matter. We can't assume that Draco has a 



steep cusp with complete impunity because many things can 
alter the inner structure of the galaxy, but under a certain 
set of conditions, such a profile becomes completely reason- 
able. In the next section, we discuss ways in which a steep 
cusp might be softened (provided that the relaxation process 
progressed far enough to form one initially), and why Draco 
should avoid this softening. 

2.3 Retaining a Steep Cusp 

Calculated as in eq. (1), i?ini„ is only on the order of a few 
AU (1 AU = 1.5 X 10^^ cm) to a few x 10^^ cm, depending 
primarily on the value of (av)^^. It stands to reason that ex- 
ternal influences such as the tidal force from the Milky Way 
should have negligible effect on a central cusp that small; 
tidal truncation is in fact observationally ruled out below a 
radius of 1 kpc (Kleyna et al., 2001b). But it is important 
that we consider any other theoretical motivation for disrupt- 
ing cusps down to such small scales as -Rmin- (The appro- 
priate scale to verify for the validity of the gamma ray limits 
derived herein is indeed i?min; however, as described in sub- 
section 3.2, the radio synchrotron emanates from within a 
larger radial size ^ 0.01 pc.) Several such ways to spoil the 
cusp are addressed here. 

(1) It has been proposed that the supermassive black holes 
found at the centers of large galaxies will change the dark 
halo cusps, although how they will change the cusps is a 
matter of some debate. One argument is that a cusp would 
steepen into a spike due to accretion of DM onto the black 
hole (Gondolo and Silk, 1999; Gondolo, 2000), causing en- 
hanced annihilation signatures. If so, then any value of the 
cusp slope index 7 over-produces the expected synchrotron 
signature beyond observational limits. Alternatively, if nu- 
clear black holes pair up in a binary system during a galaxy 
merger, as one would expect after dynamical friction pulls 
both holes to the center of the remnant but before they merge 
into one hole, then the two-body interactions of the holes 
would throw other masses out of the center (Milosavljevic et 
al., 2001; Merritt and Cruz, 2001). In that case, both bary- 
onic and dark matter would be affected, and the central cusp 
would be softened; this result is corroborated by an HST sur- 
vey (Ravindranath, Ho, and Filippenko, 2001), in that central 
mass deficits in galaxies (i.e., the departure from a cusp) cor- 
relate with the masses of their central black holes. 

The relevance of this issue depends on whether or not the 
Draco dwarf harbors a central black hole. As a general trend, 
it appears that ellipticals and spheroidals brighter than ap- 
proximately 2 X 10^° Lq do have supermassive black holes, 
and dimmer galaxies may not (Sparke and Gallagher, 2000), 
although there are exceptions. For Draco, unless its history is 
devoid of mergers, the observed velocity curve should reveal 
a softened core if there were a supermassive central hole. In 
the absence of direct evidence for such a hole in Draco, we 
assume that Draco has none, so that the concerns of the pre- 
vious paragraph only afflict other galaxies. 

(2) Could the cusp be scattered by baryonic matter? El- 
Zant, Shlosman, and Hoffman (2001) have proposed that the 
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missing cusps in galaxy cores are due to interactions between 
dark matter and clumpy baryons, in a scheme which always 
transfers energy from baryons to WIMPs, causing the bary- 
onic matter domination at small radii. Interestingly, many 
dwarf galaxies exhibit rising rotation curves out to scales like 
~10 kpc (see the discussion of dwarf galaxies in Navarro, 
Frenk, and White (1996)). But Draco, with a much higher 
M/L ratio, has constant velocity across all observed radii. 
So even if baryonic matter is to blame for destroying cusps 
in other galaxies, it has a diminished role in Draco. 

(3) Some nonlinear dynamical process among clumpy dark 
matter might expand the central region, by way of the halo 
forming from a swarm of smaller and denser subhalos which 
formed earlier. For our purposes, the effect is difficult to es- 
timate, because current simulations have not resolved much 
below Draco mass objects. However, there is no indication 
of a cutoff immediately below this mass scale in the simu- 
lations (Moore et al., 2001). Metcalf (2001) finds that even 
very small objects (> 10^ Mq, should they exist) are dense 
enough to survive tidal disruption. But these objects are too 
Ught to experience significant dynamical friction, so they still 
behave as coUisionless DM, and probably do not disrupt the 
smooth halo. Objects near the mass of the host halo will 
settle to the center due to dynamical friction (Metcalf and 
Madau, 2001), and be tidally stripped when they pass inside 
the central region (Metcalf, 2001). So is there an intermedi- 
ate scale dark object which does ruin the ~ 10^ Mq halo's 
cusp? There is no immediately obvious reason why not. 

However, it is interesting to note that globular clusters 
(typically ~ 10^ Mq) are not observed to contain signifi- 
cant amounts of dark matter (a typical globular cluster ra- 
tio of 2 Mq/Lq is given by Binney and Tremaine (1987)), 
which may mean that dark matter clumps of mass scales be- 
low those of dwarf galaxies are suppressed for some reason. 
In fact, it has been proposed that globular clusters could be 
disrupted remains of larger dwarf galaxies (Cote, West, and 
Marzke, 2001), because their numbers are appropriate to ex- 
plain the missing dwarf galaxies expected by CDM theory. 

The possibility that intermediate mass dark objects don't 
exist in large numbers can be dramatically strengthened by 
upcoming observations. If substructure DM clumps form 
on small mass scales like 10^ Mq, then thousands of such 
clumps will be observable above the cosmic microwave back- 
ground by their annihilation products (e.g., by the Planck ex- 
periment), even if they turn out to obey gently cusped NFW 
profiles (Blasi, Olinto, and Tyler, 2002). These observa- 
tions are sensitive to a lower mass cutoff on substructure 
DM clumps. So if more than ~1000 sources are not found, 
then it is difficult to see how a steep cusp could be inhib- 
ited in a dwarf galaxy like Draco by dark matter alone. (This 
holds unless the WIMPs do not annihilate like neutralinos, 
but that result is immediately consistent with non-detections 
from Draco anyway.) If these CMB foreground sources are 
found, then the results in section 4 of this paper can be con- 
strued as either limits on neutralino properties or evidence for 
important dark substructure dynamics preventing steep cusps 
in the cores of dwarf galaxies. 



3 Annihilation Signals 

3.1 Gamma rays 

Supposing that Draco has a dark SIS profile populated by an- 
nihilating neutralinos, one can calculate the gamma ray emis- 
sion caused by decaying pions made in the XX annihilation: 
tt" 77. The rate of gamma ray production above a thresh- 
old energy Eq per time per volume is 

= nl{av)^^N^{E^ > Eq) . (5) 

In this work, we adopt a method of approximation that has 
been used previously in this context (Berezinsky, Gurevich, 
and Zybin, 1992; Bergstrom, Edsjo, and Ullio, 2001) for the 
number of gamma rays produced above threshold per annihi- 
lation, Nj{E^ > Eq). The 7r° production by xX annihilation 
follows 

^ ~ K„e-^''/(x^-'^ + 0.00014) , (6) 
dx 

with X = ET^/m^c?, and constant. For the two photon 
decay process, the probability of making a photon per range 
of energy E^ is 2/£'^. Then we have 

e"^^ dx 

NJE^ > Eo) = / ^-TH 7 — , (7) 

^ ' ^JEo/m^c^ + 0.00014) X ' ' ' 

where the constant is set by requiring that one third of 
the total energy released per neutralino annihilation go into 
gamma rays, because about one third of the particles pro- 
duced by XX neutral pions, which decay via the two pho- 
ton channel. The remaining two thirds are divided equally 
among 7r+ and 7r~ particles, whose decay products will be 
the topic of the next subsection. Berezinsky, Gurevich, and 
Zybin (1992) have shown that gamma ray line flux due to 
XX ~^ 77 makes a very minor correction to the flux pre- 
dicted above, and we will ignore it here. 

To get the observable gamma ray flux (photons cm~^ s~^), 
integrate eq. (5) over the volume of the source, and divide by 
Awd^. This procedure is specific to the r-dependence of the 
halo model used, but for example, with an SIS halo where 

AA^ 1 
F^{E^ > Eo) = —^{av)^ME^ > ^o)^ , (8) 

where d is the distance from the Earth. For cusps steeper 
than 7 = 1.5, the general result is that the gamma ray flux 
varies as i?f,j7n^- 7 — 1-5^ as in the Moore profile where 
7 = 1.5, the integral over the volume of the source becomes 
dependent on an outer radius instead of Rmin- For example, 
in the Moore profile, the emission region is of size ^ and 
the flux depends only very weakly on i?min- For an isother- 
mal halo with a constant density core, eq. (8) applies, with 
^min replaced by Rcoie- Section 4 gives the results of this 
process for Draco. 
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3.2 Radio Synchrotron 



The charged pions produced in xx annihilations decay as 
tt"*" jJL^ViJt and 7r~ (jTv^ . (9) 
Muons subsequently decay via 



e^Ua^e and 



(10) 



Electrons and positrons produced in this way will generate 
synchrotron radiation if there are ambient magnetic fields. 
This calculation was first done by Berezinsky, Gurevich, and 
Zybin (1992) for the Galactic center, and is repeated in Blasi, 
Olinto, and Tyler (2002) and here, with multiple new correc- 
tions applied. (For some parts of this section, cgs units have 
been provided for clarity.) 

We need to know the number spectrum of electrons and 
positrons produced in each neutralino annihilation, dN^ / dEf, . 
We use a formulaic simplification of this function (Hill, 1983) 
which reduces to dNg/dEg ~ E^^^"^ at low energy, and 
drops toward zero as it nears the cutoff E^ 
expression used is 



m^. The full 



where r = (m^/m„ 



15 / K 



16 



-3/2 



1 - 



and 



dN, 



ml 



dE, 



E^ ml - ml 



and 



dN, 



dE,, 



E,, 



+ 



2 /E, 



3 \E,. 



(12) 



(13) 



(14) 



These last two equations, eqs. (13) and (14), give the decay 
products from charged pion and muon decays, respectively. 
Note that dN^/ dE^ without any superscripting indicates the 
number spectrum of electrons from the entire chain of decays 
following a single xX annihilation. 
The charged particle injection is 



dNe 
dEe 

-3/2 27 



(15) 



so that Qe ~ Ee'^'^r^^'' because ri;^ r~'^ for a cusp with 
logarithmic slope —7. Throughout this subsection we will 
continue to use this generic cusp. The electron distribution is 



dUe 
dEe 



= QeT, 



(16) 



where t = T{Ee,r) is the average lifetime of an electron of 
energy E^. There are four processes limiting the life of an 



electron considered in this work, and for any given combina- 
tion of E,, and r, the lifetime is the fastest of these processes. 
They are: (1) loss of energy by synchrotron radiation, (2) 
loss of energy by inverse Compton scattering (ICS) against 
the cosmic microwave background (CMB), (3) loss of en- 
ergy by ICS against the local synchrotron photons, and (4) 
annihilation with positrons. For pair annihilation, there is a 
characteristic time scale; for the other processes, there is a 
loss rate dEg/dt for which 



E,_ 



T ~ 



dE^/dt ■ 



(17) 



At a given E^, we find that e+e~ pair annihilation can 
be dominant in the center of the clump, followed outward 
by a shell dominated by ICS against synchrotron photons, 
followed by the outermost shell dominated by synchrotron 
losses and ICS against the CMB. The region dominated by 
ICS against synchrotron photons generates most of the total 
flux from the clump. Synchrotron ICS and pair annihi- 
lation are discussed below; both have the effect of removing 
emitters preferentially from the dense centred regions, reduc- 
ing the synchrotron signal. 

Employing a formula from Rybicki and Lightman (1979), 
the synchrotron power is 



dE, 



syn 



dt 



1.6 X 10 



-15 



BlEl 



erg/s , 



(18) 



for B^, = \B\ in microgauss. It is convenient to approximate 
the synchrotron as if all the radiation were emitted at the peak 
frequency (which is different for each electron): 



z/peak = 1.5 X Ki'^B^Ei Hz 



(19) 



Eqs. (15), (18), and (19) combine to give the radio signal 
we seek as follows: 



dup dEp dE,. 



erg 



dEp dv dt cm^ s Hz 



(20) 



The luminosity is the integral of ji, over the volume of the 
DM source. 

At this point, some details need to be considered. The 
above picture is correct but not complete; three amendments 
are analyzed presently. 

(Y) Inverse Compton scattering of relativistic electrons and 
positrons against the synchrotron photons created in the DM 
clump is the most important correction. ICS power is (Ry- 
bicki and Lightman, 1979) 



dEi^g 4 a2 2rr 



(21) 



where ctt = 6.65 x 10~^^ cm^ is the Thomson cross section 
appropriate in the limit hu <^ mpC^, (3 — 1, and Uph is 
the energy density in the photon field. Therefore, the photon 
population with the greatest energy density controls the rate 
of ICS. For file CMB, 



C^cmb = aT' 



cmb ' 



(22) 
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for a = 7.56 x IQ-^^ ergcm-^ R-^ and T^mb = 2.728 K; 
whereas for the synchrotron photon population, one must in- 
tegrate ji/ over all possible energies (from uieC^ to m^c^) 
and over all lines of sight. 

Unfortunately, this double integration for Ugyn requires a 
priori knowledge of the electron distribution in energy (Eg) 
and in space (r), since dUg/ dEe is contained in j^, via eq. (20). 
The electron distribution depends on the rate of ICS losses 
and depends therefore, in turn, on Ugyn (via eqs. (16), (17), 
and (21)). Consequently, our full calculation of dUg/dEe 
is numerical and iterative; however, some assumptions and 
simplifications will get us the broad brush answer as well. 

First, assume that the solution takes the form of separable 
power laws which we will solve for: 



dUg 
dK 



iEe,r)r^E° 



(23) 



Then we need a further assumption, that the dominant con- 
tribution to the photon density at some point r comes from 
the region inside r only. One might expect this to hold, since 
all photons generated inside r will pass through r, but only a 
small fraction of those produced outside r will. This is still 
an unjustified assumption which can be checked afterward 
by trying alternative values for the exponent f3; the numeri- 
cal results do validate the assumption as well. 

Let's define a shell of width Ar such that a photon spends 
Ar / c time in that region. Given these assumptions, and using 
r' as a radial variable of integration. 



1 



dE,, 



Efr'f^ 



dEsyn Ar 
dt c 



Anr''^dr' . (24) 



(The region inside i?min makes its own contribution as well, 
but it turns out to be small.) Combining with eq. (18), we get 
dUsyn/dEe ~ E^+^r>^+^. Then upon integrating over E^, 
we get Usyn ~ r' 
From eq. (17), 



./3+1 



E-'r-''- 



and from eqs. (15) and (16), 
^^(E-^/^r-'-){E-^r-^-^). 



(25) 



(26) 



So, simply equating the exponents in eqs. (23) and (26), we 
obtain a = —5/2 and /3 = —7 — 1/2, or 



dUe 

dK, 



5/2^-7-1/2 ^ 



(27) 



This holds wherever electron losses are fastest by ICS in 
the synchrotron photon field. In the case of losses being dom- 
inated by synchrotron radiation and ICS against the CMB in- 
stead, the loss mechanism is independent of r, so we would 
have diie/dEg ~ Ee^^'^r~'^^ . Since the two forms of elec- 
tron losses have different power laws in r (and the same 
power law in E,?f, there is a transition radius i?ics,syn (which 
is independent of iJg) between the two behaviors. 



For Draco, i?ics,syn — 0.01 pc, depending loosely on the 
WIMP mass and the ambient magnetic field strength. In- 
side this radius, the electron population is limited by ICS 
against synchrotron photons and decreases as r~^~^/^; out- 
side this radius the electron population is limited by syn- 
chrotron losses and ICS against the CMB, and decreases as 
r~'^^ . Interestingly, i?ics,syn becomes the characteristic ob- 
servable size of the radio source, and it is resolvable to inter- 
ferometers with baselines of order 100 km. 

(2) Pair annihilation between e~ and e+ particles would 
have an important role in determining the flux from a DM 
source in the absence of the ICS described above. But in the 
presence of ICS losses, pair annihilation makes little differ- 
ence to the resulting radio flux. 

The pair annihilation process is still important, because 
it reduces the population of the lower energy pairs, which 
would otherwise upscatter photons frequently enough to spoil 
the synchrotron spectrum. Although in much of the DM 
source, most electrons Compton scatter with photons, most 
photons do not Compton scatter with electrons because they 
escape first; it is these synchrotron photons whose survival 
we are considering now. In the case of Draco, and for the 
most optically thick line of sight through its diameter, the 
optical depth for a photon against synchrotron self Compton 
scattering is of order 10~^, making it a negligible effect on 
the resulting observable spectrum. But one needs to calculate 
the pair annihilation rate to discover that fact. 

We can estimate the lifetime for the case where elec- 
trons annihilate before they lose a significant fraction of their 
original energy via synchrotron or ICS as 

Te+ ^ 7 ^ r (28) 

with Ve — c, and the angle brackets indicating an average 
over E'e- • (This is equally vahd if we interchange all the 
plus and minus signs, but we will keep track of signs in order 
to distinguish between particles and antiparticles, for clarity.) 
With eq. (16), we have 

^ = ^el (29) 

dEe+ (ne-fTe±We-) 

The term in angle brackets is computed from the equations 
in Svennson (1982) (also consistent with Coppi and Bland- 
ford (1990)), particularly 

/•DO OC 

{av),± = /(7+) /(7-) (^)e± ^7- ^7+ , (30) 

where / is a distribution function for of Lorentz factor 7± , 
and (aiJ)e± is the angle averaged reaction rate per pair given 
explicitly in Svennson (1982); here we quote the asymptotic 
form, although the full formula has been used in the numeri- 
cal work described in this paper: 

^ (ln47+7_-2) (7+,7_»l).(31) 



((TW)e± =i 

T^eC 7+7- 

Eq. (30) is adapted by replacing f{j-) — > dng/dEf., and 
f{l+) ^ '5 (7+)' that is, for a chosen positron of energy 7+, 
its survival time is a function of the electron distribution. 
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Then, we have {Uf.- (je± Ve-) expressed as an energy inte- 
gral which depends on the distribution that we're trying 
to solve for. Therefore this pair annihilation calculation must 
be done numerically. The dependences here are nested in- 
finitely, which is to say that Tg+ depends on n^- which de- 
pends on Te- which depends on ne+ and so on. It is most 
efficient to try to determine the correct dependences a priori, 
and then integrate accordingly. 

So if we neglect the natural logarithm factor built into 
((T?J),.± from eq. (31), because it evolves slowly compared 
with powers of 7+ and 7_, then we can proceed to find the 
energy dependence for Tg+ to within a factor of about E^°'°^ 
in accuracy. 

The reaction rate for a given positron is (ng-o-g±Wg-), 
which varies inversely with electron energy (any other hy- 
pothesis leads to a contradiction). Therefore, the positron 
will preferentially annihilate against an electron at the low 
end of the energy distribution, so that in eq. (29), ((Te± Ug- ) ~ 
E~^, and rig- is roughly constant. That is, a given positron 
of any energy sees essentially the same distribution of tar- 
get electrons, so to a first approximation, it is sufficient to 
replace the electron distribution with a single population at 
energy E^- ^nieC?. By that rationale, Tg+ ^ Eg+. 

With this last result in hand, it becomes useful to use 



per-unit-length absorption coefficient: 



Tg+ 



Ee- 

Ep+ 



(32) 



when doing the full numerical integral over target electron 
energies, to avoid the infinite nesting problem. At low ener- 
gies where dNg/dEe ~ Eg^^'^, this means that dng/dEg 
E~^/^ and rig ^ E^^^ (where the subscript e refers to either 
an electron or a positron). 

The distribution in space, rather than energy, is easier to 
determine. In order to satisfy eq. (29) with the injection qg ~ 
r"^''' set by a cuspy DM profile, we must have dUg/dEg ~ 
r~'^ . (As before, —7 is the halo profile logarithmic slope, 
while 7+ and 7_ are Lorentz factors.) Then for a region 
whose losses are fastest by pair annihilation, we have 



dEg 



(33) 



So now we have another transition radius between the r 
pair annihilation distribution and the r^^~'^/'^ ICS distribu- 



tion, which we designate -Rann,ics- If -Ra 



> iJmin, then 



we add a pair annihilation dominated central region to the 
picture described above regarding ICS. But for most choices 
of and Draco's magnetic field strength, iiann.ics < Rmin 
for electrons radiating at frequencies of interest, and the dis- 
tribution flattens at the center before the density can get high 
enough to create a pair annihilation dominated region. 

(3) Synchrotron self absorption (SSA) is where an elec- 
tron in a magnetic field can absorb the synchrotron photon's 
energy. Still linking photon frequency to the energy of its 
emitting electron with eq. (19) for simplicity, the appropri- 
ate adaptation from Rybicki and Lightman (1979) yields the 



dE, 



a,. 



syn 



Stti/ dt 



d_ 
dv 



1 dug 
V dv 



(34) 



SSA absorbers of a particular photon are electrons of similar 
energy to the original synchrotron emitter which made the 
photon. As such, SSA is only effective below some i/grit. 
where the corresponding electrons are more numerous. 

The source function for an absorbing source is = jv/ ctv, 
and the optical depth is 



'^v = j OL^, 



dz , 



(35) 



where 2: is a Une-of-sight coordinate. The flux density from 
such a source at distance d from the Earth is obtained by 
summing over each line of sight through the DM clump: 



erg 



cm^ s Hz 



(36) 



The SSA cutoff frequency Vgrit depends in part on the ab- 
sorber density, and without ICS or pair annihilation ef- 
fects, SSA would be the most important correction to the 
calculated flux in a steep cusp. SSA has been included in 
the calculations performed here, but since the central ab- 
sorber population is diminished by ICS, it is of small im- 
portance. Using appropriate parameters for Draco, we find 
fcrit — 10 MHz, whereas the limiting observations discussed 
in the next section are at 4.9 GHz. (It is fortunate that SSA 
has little impact on the problem, because otherwise it would 
create serious complications in the calculation of ICS against 
the synchrotron photon field.) 

4 Results 

Supersymmetric dark matter has a number of variable param- 
eters which affect its annihilation and clustering properties. 
However, these parameters boil down to only two relevant 
quantities for annihilation signature searches, and those are 
and (av)^^. Given external constraints on a dark mat- 
ter source (such as Draco's total mass and the matter density 
local to it, in the present case), the mass determines the 
X number density. The more massive the particle, the fewer 
are needed, so the annihilations are less frequent. We present 
here the collection of and {cfv)-^^ values available to the 
neutraUno, derived from observations of Draco. 

Before proceeding to these limits, we note for compar- 
ison another indirect detection approach considered in the 
literature, that of detecting neutrinos (as from eqs. (9) and 
(10)) coming from neutralino annihilations at the center of 
the Earth and the Sun. By elastic scattering with nuclei in the 
Sun or Earth, a WIMP can lose enough energy to be gravi- 
tationally captured by that body. This captured population 
is the source of the neutrino signal in question. Supersym- 
metric parameter space can be ruled out in this way by neu- 
trino telescopes such as Baksan or AMANDA (for example, 
see Bergstrom, Edsjo, and Gondolo (1997) and Ahrens et al. 
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(2002)). This method can currently eUminate more parame- 
ter space than direct detection experiments, and less than the 
Draco EGRET limits as calculated herein. The constraints 
from Draco require the extrapolation of its halo density pro- 
file to small radii; the constraints from neutrino telescopes 
instead require the additional steps associated with elastic 
scattering, capture, and neutrino propagation through mat- 
ter. Of course, this approach and the one used in the present 
work are both indirect detection scenarios; for a review of 
direct experimental searches, see e.g. Morales (2001). 

4. 1 Gamma Ray Constraints 

Gamma ray data specifically on Draco are lacking, so the 
current limiting observations come from an all-sky survey. 
The EGRET (Energetic Gamma Ray Experiment Telescope) 
instrument, flying on the CGRO satellite (Compton Gamma 
Ray Observatory), is a pair production telescope (see Mattox 
et al. (1996) for an introduction to EGRET and its data analy- 
sis). As such, its useful beam width is around 30°, making it 
well suited for a full sky survey; its energy range is ~30 MeV 
to below ~100 GeV (Sreekumar et al., 1998). 

EGRET sources are usually counted by their photons above 
100 MeV; however. Lamb and Macomb (1997) have per- 
formed a separate binning of the data, making a catalog of 
gamma ray sources above 1 GeV, which is of particular inter- 
est here. The least significant detection (4cr) in the catalog, 
with the least flux, is the Large Magellanic Cloud (LMC), 
at (1.1 ± 0.4) X 10"^ photons cm"^ s^^ above 1 GeV. We 
therefore take 10~* photons cm~^ as the flux limit for 
detection, and require that Draco emit less. 

In fact, since this limit is only based on non-inclusion 
in the GeV catalog, specific knowledge of the exposure on 
Draco could result in a bound somewhat better that the one 
derived herein. A subsequent analysis of the EGRET data 
provides an estimate of the instrument's extragalactic gamma 
ray background (Sreekumar et al., 1998). With Draco located 
well outside the Galactic disk (I = 86.37°, h = 34.71° (van 
den Bergh, 2000)), this extragalactic component is the pri- 
mary background source. The extragalactic flux, adapted to 
the GeV catalog's 30 arcmin pixel size, 



F-^{E-^ > Eo) = 8.6 X 10" 



Eq 

GeV 



photons 



(37) 



is well below the level of one photon per minimum EGRET 
exposure. We therefore separately consider this background 
in our calculations, conservatively assuming that EGRET's 
exposure on Draco is their minimum exposure (roughly true), 
for the case of no photons detected within 30 arcmin. We 
present this as an example of the type of improvement possi- 
ble without coUecting new data. 

In figure 1, we plot the neutralino parameter space. In 
this figure, we have presented limits derived by use of the 
SIS halo profile with different curves giving different present 
and future observational bounds; a comparison of halo mod- 
els will follow in subsection 4.3. Both EGRET constraints 




allowed region 



EGKEI' (E„=l GeV) 

EGRET (with eg bg) 

GLAST (E, = l GeV) 

— VERITAS (E| = l()() GeV) 
VERITAS (E =1 TeV) 



10 10 10 

in GeV 

Fig. 1. The — {(t«)xx plane for neutralino dark matter. Curves 
indicate results where the SIS halo profile has been used. Shown 
as solid lines are the gamma ray constraints from non-inclusion in 
the EGRET all-sky survey, where the thick line is the current bound 
and the thin line is a potential improvement using an estimate of 
the extragalactic gamma ray background (see text). Both employ 
a threshold Eo — 1 GeV. Also shown are two dashed lines show- 
ing upcoming VERITAS constraints with cutoffs at 100 GeV and 
1 TeV, and a dotted line for GLAST at Eo = 1 GeV (sensitivities 
for all of these taken from Weekes et al. (2001)). The region below 
the curves is allowed by gamma ray observations. Uncertainties in 
Draco's mass and distance lead to at most a factor of 3 change in 
the vertical axis position of these curves. The dots denote typical 
supersymmetric dark matter properties (Gondolo et al., 2000). 



are shown: the current bound (thick solid Une), and the ex- 
ample bound with the newer extragalactic background (thin 
solid line). Points above the curves are prohibited by these 
observational constraints from Draco. 

Two notes should be made regarding accuracy. First, the 
curves in the figure assume that MDraco = 8.6 x 10*^ Mq, 
which is measured using observations of its luminous mat- 
ter (Kleyna et al., 2001a). So this mass samples the inner 
halo, and the total mass could be higher. In this sense, the 
curves plotted should be considered conservative. Second, 
recall from the discussion of i?inin in subsection 2.1 that a 
higher estimate of i?min would result from a slower relax- 
ation process. At most, this would make these curves less 
constraining (i.e., raise them) by a factor of about 5. 

The dots in the figure survey the possible combinations of 
and {<tv)^^ for neutralino dark matter (Gondolo et al., 
2000). The dotted region is characteristic of supersymmetric 
parameters often considered in the literature (e.g., Gondolo 
and Silk (1999); Baltz et al. (1999)). They are constrained by 
the requirement that the relic density be cosmologically inter- 
esting; that is, large enough to account for a significant com- 
ponent of the matter density without exceeding it. The choice 
of these "interesting" values is somewhat arbitrary and au- 
thors differ slightly, but a restriction like 0.05 < ft^h"^ < 0.5 
as depicted in the figure is typical. It should be noted that 
values of < 38 GeV are disfavored by accelerator exper- 
iments (Serin, 2000). 
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Fig. 2. The — {ov)^^ plane for neutralino dark matter. Curves 
indicate results where the SIS halo profile has been used. Shown as 
thin lines are the radio synchrotron constraints from VLA observa- 
tions. From top to bottom, these lines correspond to magnetic fields 
of 0.01, 0.1, 1.0, and 10.0 in the core of Draco. The region 
below these curves is allowed by radio observations. Uncertainties 
in Draco's mass and distance lead to at most a factor of 3 change 
in the vertical axis position of these curves. The dots denote typical 
supersymmetric dark matter properties. 



Future gamma ray telescopes offer better constraints. The 
upcoming satellite mission GLAST (Gamma ray Large Area 
Space Telescope), will offer much improved constraints as 
shown in the figure (Gehrels and Michelson, 1999). It is ex- 
pected to launch in 2005 for a two year all-sky survey, effec- 
tive between 20 MeV and 300 GeV, with a sensitivity about 
20 times better than EGRET at 1 GeV. As can be seen in 
the figure, GLAST has the potential to rule out most of the 
neutralino DM parameter space. 

Apart from satellite missions, atmospheric Cherenkov tele- 
scopes (ACTs) are ideal for measuring at high energy. Rather 
than detecting the photons via interactions inside the tele- 
scope, ACTs count gamma rays by the air showers they pro- 
duce in the Earth's atmosphere. Specifically, they respond to 
Cherenkov Ught produced by charged particles in the shower 
traveUng faster than the speed of light in air. The next gener- 
ation gamma ray telescope VERITAS (Weekes et al., 2001) 
can provide significant benefits for this study with a typical 
exposure on Draco (we quote here an exposure of 50 hours). 
Curves for VERITAS with threshold energies of 100 GeV 
(long dashes) and 1 TeV (short dashes) are depicted in the 
figure. VERITAS is sensitive between about 50 GeV and 
50 TeV; it is expected to see first Ught in 2002 and become 
fully operational in 2005. 

4.2 Radio Constraints 

The applicable radio continuum limits for Draco were per- 
formed by Fomalont and Geldzahler (1979), using the VLA 
(Very Large Array) facility in New Mexico. At 4.9 GHz, they 
report no flux from Draco above 2 mJy, at the 3(7 level. Their 
search spanned a radius of 4 arcmin around the galaxy center. 



Two additional assumptions beyond those discussed in sec- 
tion 2 are needed in order to place a radio limit, because syn- 
chrotron emission requires magnetic fields. 

First, we must assume a magnetic field strength \B\ for 
the core of Draco. No such estimates exist in the literature, 
so this paper must be viewed as constraining either the neu- 
tralino or Draco's magnetic properties. To give some bearing 
on the estimation of we can consider measurements for 
other dwarf galaxies. The Magellanic Clouds, for example, 
each carry fields of ~5 /zG (Pohl, 1993). A survey of low sur- 
face brightness dwarf galaxies (similar to Draco) leads Klein 
et al. (1992) to infer typical field strengths between 2 and 
4 fjbG, although this survey was conducted near 5 GHz and 
Draco specifically gives no measured signal there. The dwarf 
irregular NGC 4449 was found to harbor ~14 /zG fields de- 
spite the lack of ordered rotation (Chyzy et al., 2000), sug- 
gesting that some non-dynamo process may be sufficient to 
generate microgauss fields. Based on these data, it seems 
appropriate to expect ~ 0.1 — 1 fiG in Draco, although at 
present there is no compelling evidence. 

Second, not only must there exist magnetic fields, but they 
also must be sufficiently constricted, so as to trap electrons 
and positrons long enough for them to radiate. Implicit in 
the computation of flux from DM sources in this paper is the 
assumption that particles never stray far from their birth 
places — an assumption of strong magnetic confinement. To 
begin contemplating this, consider that the gyroradius of any 
electron relevant to this problem is significantly smaller than 
the smallest length scale applicable; that is, R^yro ^ ^min- 
So if we assume that electrons are reflected back and forth 
along a field line every 10 — 100 i?gyro> then strong con- 
finement is justified. Klein et al. (1991) notes that mag- 
netic fields do in fact appear to be highly disordered in some 
dwarf galaxies (Draco not considered), based on lack of ra- 
dio polarization observed. The confinement time needs to be 
longer than the inverse Compton time scale against the syn- 
chrotron photon field in order to validate our model. This 
duration is typically around 10^ — 10^ years for key emitting 
regions in Draco's halo, depending on neutrahno parameters, 
for 4.9 GHz radiation as Umited by the VLA. 

These caveats noted, the synchrotron constraints are plot- 
ted in figure 2, where again the points above each curve can 
be ruled out by it. Limits for 0.01, 0.1, 1, and 10 /uG fields 
are shown here for the SIS halo; this choice will be varied in 
section 4.3. As with figure 1, the possibility that Moraco was 
underestimated makes these curves conservative. 

Although the radio limits are less direct because of as- 
sumptions that must be made about the B field, they can be 
stronger than current gamma ray limits. For example, if the 
field is found to be at least 1 /iG and adequately confining, 
then (a-v)xx for ^ neutralino is restricted to be be- 

low 4 X 10~^^ cw?/s, as compared with ~ 10~^* cm^/s 
in gamma rays. Heavier neutralinos are strongly constrained 
at 0.1 jiG, but they aren't limited at all by current gamma 
measurements. On the other hand, if the magnetic fields are 
weak, then synchrotron offers no constraint on light x parti- 
cles. For example, a 100 GeV neutralino's annihilation prod- 
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Fig. 3. Parameter space constraints from the EGRET 1 GeV catalog 
with seven different halo models for Draco. The halos are: cusps 
with 7 = 2.0 (solid line), 1.9 (dot-dash), 1.8 (dot-dot-dash), and 1.7 
(dot-dash-dash), Moore (dotted line), SIS with constant density core 
at 1 pc (thick dashed line), and SIS with constant density core at 
0. 1 pc (thin dashed line). The thickest solid line is repeated from the 
current EGRET bound given in figure 1. Lines not shown indicate 
that the corresponding halo yields no constraints on this plot. 

ucts aren't energetic enough to radiate at 4.9 GHz (where the 
VLA Umit is) if the field is ^^0. 1 /xG or below. 

As in figure 1, the dots indicate the SUSY parameter space. 
Here we see the importance of the magnetic field. At 1 /j,G 
one can exclude the majority of neutralino dark matter can- 
didates; at 0.01 /iG one can exclude almost none of it. 

One further point worth noting is that the gamma ray and 
radio analyses, taken together, permit us to constrain the mag- 
netic field properties of Draco. If a future gamma ray detec- 
tion is made, then the coiTesponding line in the — (trw)^^ 
plane can be transposed to figure 2, with the result that some 
values of will not be allowed. 

4.3 Other Halo Profiles 

Figures 1 and 2 solely employ the SIS to generate the Umiting 
curves. In this subsection, we consider the sensitivity of our 
results to the choice of halo profile. 

In figure 3, we present modified gamma ray results cor- 
responding to the current EGRET bound shown as the thick 
line in figure 1 . All curves on this plot correspond to the same 
EGRET observations above the Eq ^ \ GeV threshold, but 
for different halo models. Figure 4 presents the same halos, 
but for the GLAST curve from figure 1 with Eq — \ GeV. 

The halo models considered all assume spherical symme- 
try. They are: cusped profiles with 7 = 2.0 (SIS), 1.9, 1.8, 
and 1.7, the Moore profile as described in subsection 2.1, 
and the SIS truncated with constant density cores of radius 
1 pc and 0.1 pc (the velocity dispersion data extend down 
to ~ 10 pc). Figures which omit some of these curves do 
so because those curves do not constrain any of the param- 
eter space shown. For example, with Draco modeled by the 
Moore profile, the EGRET data are insufficient to contrain 



Fig. 4. Same as figure 3, but now showing modifications to the 
GLAST line from figure 1, rather than EGRET. 

{(jv)^^ below ^ 10~^^ cm^/s; one can do somewhat better 
in the future with GLAST (figure 4). 

For the radio case, the corresponding results are presented 
in figure 5 for \B\ — 0.1 /iG, and in figure 6 for \B\ = 1 /iG. 
In each figure, the same seven halo models are shown as in 
figures 3 and 4, including the appropriate SIS curve from 
figure 2. 

The cusps 7 = 2.0, 1.9, 1.8, and 1.7 are calculated with the 
full inverse Compton scattering coiTections described in sub- 
section 3.2. Applying eq. (27) yields dn^/dEf. dependences 
which go as r^'^-^, r^^ *, r~'^-^, and r~^-^, respectively. 

Inverse Compton scattering does not play a significant role 
in the SIS core cases or in the Moore case, because in each of 
these halos, the primary emitting region is outside of i?ics,syn- 
The Moore case turns out to create a diffuse enough collec- 
tion of electrons and positrons that synchrotron self absorp- 
tion is also negligible, although SSA does significantly affect 
the results for the isothermal core models. 

It is interesting to note that for both gamma ray and radio 
constraints, the Moore profile yields more flux than a j — 
1.7 cusp, even though the inner Moore halo only has 7 = 
1.5. This is a result of the fact that the Moore profile does 
not extend as r~^-^ indefinitely, but rather is concentrated 
into Tg. This concentration effect makes up for the lower 7 
exponent. 

From figures 3 through 6, it is clear that the choice of halo 
model is of great importance in determining the resulting 
dark matter constraints. Even p ~ t-^I-S yields weak lim- 
its with EGRET data. With GLAST, many of the steeper 
models are valuable, but less dense ones are not. In the ra- 
dio, a magnetic field of 0.1 pG is roughly as constraining 
at high values as GLAST is at lower values. At 
1 /iG, the situation is much better, with all examined pro- 
files generating strong constraints. In this case, the Moore 
halo too constrains WIMPs well, although the Moore profile 
is not particularly constraining in any of the other cases con- 
sidered. The primary lesson from these example cases is that 
the results do vary appreciably even among relatively steep 
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Fig. 5. Parameter space constraints from the VLA with seven dif- 
ferent halo models for Draco, assuming a magnetic field strength 
\B\ = 0.1 fiG. The halo modifications are the same as those in 
figures 3 and 4. 

cusps, underscoring the importance of pinning down Draco's 
central density structure accurately. 

5 Summary 

This paper applies gamma ray and radio observational limits 
to the behavior of dark matter in the Draco dwarf galaxy, only 
79 kpc distant from the Earth. Recent measurements of stel- 
lar velocities in Draco suggest that it is extremely dark mat- 
ter dominated, and that its halo distribution is nearly isother- 
mal. Because such a halo should be densely concentrated at 
the center, we expect WIMP annihilation to proceed rapidly 
there, producing observable gamma ray and synchrotron ra- 
diation. 

Current limits derived herein impose restrictions on the na- 
ture of the dark matter particle. If it is a neutralino (or simply 
annihilates into pions like one) then figures 1 through 6 con- 
strain the mass and annihilation cross section available to it. 
These results emphasize the importance of the choice of halo 
model for Draco. Our strongest constraints emerge with the 
SIS profile. In this case, EGRET data restrict much of the 
supersymmetric parameter space, and if Draco has a suitable 
magnetic field, then the VLA data can exclude significantly 
more of it. Other halo profiles with a core or a less steep 
cusp generally become significantly constraining with future 
gamma ray experiments, or in the radio with magnetic fields 
in the fiG range. 

An important conclusion from this work is that further par- 
ticle dark matter research would be well served by increas- 
ingly detailed studies of Draco. 
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Fig. 6. Same as figure 5, but now showing modifications to the 
= 1 line from figure 2. 
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